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R. C. Whitten, W. J. Borucki, I. G. Poppoff,* L. Glatt,t G. F. Widhopf,t L. A. Capone,t and C. A. ReigelJ 

Ames Research Center 


An earlier assessment (NASA Reference Publication 1026) of potential effects of advanced supersonic aircraft 
on stratospheric ozone employed one-dimensional photochemical models to obtain estimates of globally averaged 
changes in ozone. The study reported here, an extension of the earlier work, was undertaken to estimate the lati- 
tudinal dependence of the ozone perturbation; two two-dimensional photochemical models were used for the 
purpose. Two SST design concepts assumed in the previous work were also employed for the work reported here. 

As predicted in an earlier study, it appears that realistic SST fleet sizes should not cause concern with respect 
to the reduction of the total ozone overburden. For a fleet of 250 aircraft, the change in the ozone column is 
predicted to be very close to zero; in fact, the ozone overburden may actually increase as a result of nitrogen 
oxides deposited at altitudes of 1 7.5 or 20 km. The calculations performed, using the chemical system described 
herein (ca, 1978), show that above 25 to 30 km the ozone abundance decreases via catalytic destruction, but at 
lower heights it increases, mainly as a result of coupling with odd-hydrogen species. Water vapor released in the 
engine exhaust is predicted to cause ozone decreases; for the hypothetical engines used in the study, the total 
column ozone changes due to water-vapor emission largely offset the predicted ozone increases due to NO^ 
emission. The actual effect of water vapor may be less than calculated because present models do not include 
thermal feedback. ^Feedback*' refers to the cooling effect of additional water vapor that would tend to slow the 
HOx reactions which destroy ozone. 

Because it is also appropriate to the study reported here, the reader is referred to the last paragraph of the 
^"Summary ” of the earlier work. 


I. INTRODUCTION 

In a previous publication (Poppoff et al., 1978, 
hereinafter referred to as “Part I”) the results of an 
assessment of the potential effect of future aircraft 
operations on the abundance and distribution of 
ozone in the stratosphere was reported. The assess- 
ment was made with a one-dimensional photo- 
chemical model (Turco and Whitten, 1977). It was 


found that in contrast to the earlier predictions in 
the Report of Findings of the DOT’s Climatic Impact 
Assessment Program (Grobecker et al., 1974) and the 
National Academy of Science report, “Environmental 
Impact of Stratospheric Fli^t” (1975), nitrogen 
oxide emissions by aircraft cruising in the strato- 
sphere could, for some cases of injection, cause 
increases in the ozone column density. i The situation 
appears to be reversed with respect to water vapor 


♦Presently at Flick Point Associates, Carnelian Bay, Cali- 
fornia. 

tAerospace Corporation, El Segundo, California. 

JSan Jose State University, San Jose, California. 


1 Ozone depletions aie still found at higher altitudes 
in the stratosphere; however, ozone increases in the lower 
stratosphere are greater than the upper stratosphere deple- 
tions. 




emissions by aircraft engines: the older work pre- 
dicted a negligible effect of water vapor on strato- 
spheric ozone whereas the recent assessment in Part I 
led to predictions of small ozone losses. Still, the 
assessment of the water-vapor effect is in some doubt 
because of the neglect of thermal-photochemical 
feedback (e.g., Luther and Duewer, 1977), which 
tends to diminish the water-vapor effect. Injection 
of water vapor in the stratosphere may also cause a 
“greenhouse” increase in surface temperature and a 
decrease in stratospheric temperature. Based on the 
computations of Manabe and Wetherald (1967), we 
estimate these temperature changes to be ^^0.04° C 
and 1.0*^ C, respectively. 

These changes in our ideas about potential effects 
on the ozone layer of aircraft cruising in the strato- 
sphere are due largely to relatively recent direct 
measurements (Howard and Evenson, 1977; Burrows 
et al., 1977) of the rate coefficient for the reaction 

NO + HO2 ^ NO2+OH (1) 

The newly measured value is about 3040 times 
larger than the older one obtained by indirect means 
and has altered many of our ideas of stratospheric 
aeronomy. The aeronomical effects of a large rate 
for reaction (1), discussed in detail in Part I, are 
summarized later in this report. It is worth mention- 
ing at this point that a recent direct measurement 
(Zahniser and Howard, 1978) of the rate coefficient 
of the reaction 

HO2 + O3 ^ OH + 2O2 (2) 

leads to predictions of even larger ozone increases due 
to NOx emissions. The work of Zahniser and Howard 
(reaction (2)) was reported too late to be used in 
Part I and for most of the computations reported 
here; however, some computations were made with 
the new value of the rate coefficient and are 
presented for comparison. 

Assessments using two- and three-dimensional 
models of the latitudinal variation of the effect of 
aircraft operations on ozone have been reported 
(Alyea et al., 1975; Borucki et al., 1976; Cunnold 
et al., 1977; Widhopf et al., 1977; and Hidalgo and 
Crutzen, 1977). They all used values of the rate 
coefficient for reaction (1) that were considerably 
smaller than those of Howard and Evenson (1977) 
and Burrows et al. (1977). As would be expected, 
all of these calculations predicted ozone losses 
due to supersonic transport operations, although 


Widhopf et al. (1977) and Hidalgo and Crutzen 
(1977) did predict that ozone increases would be 
caused by aircraft operating below about 14 km. 
Using the new value of the rate coefficient for reac- 
tion (1), Widhopf and Glatt (1979a) and Glatt and 
Widhopf (1978) predicted ozone increases in the 
stratosphere due to combined NOx and H 2 O emis- 
sion. It should be noted that these differences in 
the results reported from the earlier studies and 
those reported here are due almost entirely to the 
use of the earlier, smaller values for the rate coeffi- 
cient (1) and not to any basic lack in the models or 
the theory used by the groups making the assess- 
ments; however, there have been refinements in 
both theory and models since the earlier studies. 
In section V we will return to the results of the 
earlier studies. 

A new study of potential aircraft effects on the 
ozone layer, using multidimensional models, is thus 
needed in order to bring the assessments up to date. 
It is the purpose of this report to try to meet this 
requirement with the aid of two two-dimensional 
models, one developed at Ames Research Center 
(Borucki et al., 1976; Whitten et al., 1977) and the 
other developed at the Aerospace Corporation 
(Widhopf etal., 1977).^ The models used the same 
reaction rate coefficients and similar photolytic 
cross sections, but different numerical techniques 
and different transport parameterizations; although 
the new value of the rate coefficient for reaction (2) 
was reported too late for incorporation into the 
Aerospace calculations for this particular study, we 
have been able to include it for a few of the Ames 
computations. A detailed comparison of the model 
results and a determination of the specific reasons 
for model differences leads to added confidence in 
the results of both models. 

In the following sections, we discuss our model 
of aircraft emissions (sec. II), the two photochemical 
models used in the study (sec. Ill), predictions of 
trace constituent ambient distributions (sec. IV), 
and the results of the assessments (for which the 
calculations were performed mainly in 1978 and 
completed in January 1979), a comparison with 
earlier work, and discussion of the assessments in 
terms of stratospheric aeronomy (sec. V). Finally, our 
results are summarized in section VI. 


2 Developed in part under funding from the FA A High 
Altitude Pollution Program, U. S. Department of Trans- 
portation. 
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IL EMISSIONS MODEL 


In section II of Part I (Poppoff et al., 1978) we 
discussed the potential emissions of nitrogen oxides 
from two suggested SST design concepts (denoted 
types A and B), Type A, documented in a report by 
Baber and Swanson (1976), is based on the AST-100 
airframe and on a single-spool turbojet with variable 
turbine geometry. The flight conditions assumed in 
the document include a speed of Mach 2.7 and a 
cruise profile varying from 19 to 21 km, nominally 
centered at an altitude of 20 km. The initial fuel con- 
sumption rate for this particular case is 44,100 kg/hr 
and the average consumption rate over the cruise 
portion of the flight is 37,800 kg/hr. For the type B 
aircraft, a set of conditions that are representative 
of design concepts using the Pratt and Whitney VSCE 
502B variable-cycle study engine and more advanced 
airframes was used. For these conditions, the cruise 
speed is Mach 2,3, the cruise profile is centered at an 
altitude of 17.5 km, the initial fuel consumption rate 
is 41,050 kg/hr, and the average fuel consumption 
rate is estimated to be 35,200 kg/hr. The essential 
difference between the two design concepts (insofar 
as potential effect on stratospheric constituents is 
concerned) is the flight altitude. In the studies 
reported here, far-term engine emission technology 
for both type A and type B aircraft was assumed 
(6 and 7 g N 02 /kg fuel, respectively; see table 1 
of Part I). Water vapor is also emitted by aircraft 
engines; the amount is determined stoichiometrically 
in well-burned fuel and is equivalent to an emission 
index of about 1.3 kg HaO/kg fuel. 

In modeling the latitudinal distribution of aircraft 
emissions, zonal symmetry was assumed, that is, that 
the effluents undergo zonal spreading rapidly and 
uniformly. This approximation is a reasonable one 
because of the strong zonal wind fields (Belmont 
et al., 1975). A hypothetical route structure was 
then constructed by selecting from the April 1977 
Official Airline Guide (Worldwide edition) those 
long-range (^2000 n. mi.) over- water routes that 
were then (1977) flown by Boeing 747’s (including 
747-SP). Such B-747 traffic was selected because it 
was believed to be representative of the type of 
traffic that might be partially transferred to large 
(250-300-passenger) supersonic transports in the 
future. In constructing our traffic model the mag- 
nitude of future Southern Hemisphere operations 
may have been underestimated, since they may grow 


faster than the Northern Hemisphere counterpart; 
unfortunately, there are no quantitative estimates of 
such asymmetric growth. 

Having established a route structure, the appro- 
priate great circle routes and the emissions in each 
latitude increment (5^^ for the Ames model, 10"^ 
for the Aerospace model) were computed. The 
aircraft were assumed to rise in altitude linearly 
as fuel was consumed, with nominal midcruise 
altitudes of 20 and 17.5 km for types A and B 
aircraft, respectively (i.e., the altitude spread for a 
4000-n. mi. flight by type A aircraft is 19 to 21 km 
and 16 to 18 km for the type B aircraft). Fuel burned 
in reaching cruise altitude was not included because 
the “flight envelope” of the consumption, that is, 
the rate of consumption as a function of altitude, 
was not known. From the design data reported by 
Baber and Swanson (1976), the amount of fuel so 
consumed in the stratosphere appears to be rather 
small compared with fuel burned at cruise altitude. 
The NOx emission rates for each altitude and latitude 
increment were then adjusted such that the simulated 
global emission rate was that characteristic of 250 
type A or type B aircraft flying 7 hr/day in the cruise 
mode. The NOx emission data for type A aircraft 
with an NOx emission index of 6 g N 02 /kg fuel are 
shown in figure 1. The same distribution was used 
for the type B aircraft simulation by lowering the 
altitude regime accordingly. 

The effect on stratospheric ozone of the combined 
engine emissions of supersonic (type A far-term 
technology) and subsonic aircraft is also considered 
in this report. For the subsonic aircraft emission 
model, the NOx emission rates suggested by Oliver 
et al. (1977) were used; the subsonic emission esti- 
mates are presented in table 1. 

III. THE PHOTOCHEMICAL MODELS 


It is difficult to model the atmosphere, especially 
the troposphere. Atmospheric models necessarily use 
a finite grid spacing (or a finite number of waves in 
spectral representation) and thus ignore the coupling 
of sub grid scale motions with the large-scale motions 
that can be modeled. Two-dimensional models are 
more restricted than three-dimensional models 
because they can simulate only zonally -averaged 
atmospheric conditions and thus cannot follow the 
development of planetary wave systems. Instead, it 
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is customary to specify empirically a set of transport 
parameters composed of those that specify a large- 
scale meridional bulk motion (usually by a stream 
function) and those that specify smaller-scale pro- 
cesses that lead to mixing (by “eddy transport” 
coefficients). The bulk motion can be calculated by 
zonally averaging the wind fields predicted by three- 
dimensional models or by computing the meridional 
wind field from the energy and momentum equa- 
tions (e.g., Vupputuri, 1975; Harwood and Pyle, 
1977). Eddy coefficients are then obt^ed indirectly 
from measurements and adjusted if necessary until 
the computed distributions match the observed 
distributions of atmospheric tracers, such as those 
discussed in section IV. Two-dimensional models 
are thus highly empirical in nature, and they are 
incapable of simulating the observed variability of 
atmospheric constituents. 

Because the two photochemical models employed 
for the current study have been described in detail 
in the literature (Borucki et al., 1976; Whitten et al., 
1977; Widhopf et al., 1977), they are only briefly 
described here. Both models are based on the govern- 
ing equations for trace constituents in a meridional 
plane 





( 3 ) 


where P- and are production and loss rates, respec- 
tively, for the /th constituent, v is meridional bulk 
velocity, <I> is the “eddy flux” which represents 
large-scale eddy motions, and is the number 
density of the /th constituent where the Ames model 
solves directly for n- and the Aerospace model com- 
putes the corresponding m^ng ratio; and L- 
will, of course, have different forms, depending on 
whether equation (3) is written in terms of the 
number density or mixing ratio. For example, the 
“eddy” fluxes in each model are expressed as 




-K^ " model) 

[vW] 


( 4 ) 

(Aerospace model) 


in which represents the eddy-diffusivity tensor 
(vertical, meridional, and “off-diagonal” terms); 
p, H, and T are the atmospheric density, scale height, 
and temperature, respectively; and is a vertical 
unit vector. 


Ames Model 

Spherical geometry is employed in equations (3) 
and (4) with the horizontal coordinate extending 
along a meridian from lat. 80° S to lat. 80° N in 5° 
intervals and the vertical coordinate extending 
from the ground to an altitude of 60 km in 2.5-km 
intervals. The end boundary conditions assume zero 
flux of all constituents across the vertical boundaries. 
The upper boundary conditions are specified by zero 
vertical flux for all species, and the lower boundary 
conditions are either (1) chemical equilibrium or 
(2) fixed concentration that is usually representative 
of measured values. Species having very rapid loss 
mechanisms (e.g., O, OH, HO2) were assumed to 
be in photochemical equilibrium at the lower 
boundary. Long-lived species were held fixed at 
the ground (e.g., N2O at 7.5X10^^ cm”^, O3 at 
4.5X10^^ cm”^). For the computations reported 
here, we fixed the concentrations of the water soluble 
species HNO3, HCl, and H2O2 equal to zero at the 
lower boundary. The mixing fractions of CO2 were 
fixed at all grid points at 367 ppmv. 

The chemical rate equations are solved using an 
implicit technique. Omitting the transport term from 
equation (3) yields the finite difference form (with 
Q-P-L) 




A/ 




( 5 ) 


which can be linearized by taking the first term in a 
Taylor series expansion about Q- : 


/+! 



A/ 




k 


( 6 ) 


The members of the set of equations are coupled and 
require solution by matrix methods; the “families” 
approach was not used for the computation reported 
here. Diurnal averaging of photodissociation rates 
and reaction rate coefficients is performed, using, 
respectively, the techniques reported by Cogley and 
Borucki (1976) and Turco and Whitten (1978). The 
transport computations are time-split into those for 
vertical and horizontal advection and diffusion. 
Solutions for the transport processes are obtained 
from a mass-conserving, forward-time, space-centered 
finite difference formulation of the equations. An 
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implicit formulation is used for the horizontal diffu- 
sion, and an explicit formulation is used for all other 
transport processes. 

The atmospheric meridional bulk velocity and 
eddy transport fields are derived in the manner 
discussed by Borucki et al. (1976) and Whitten 
et al. (1977). The eddy coefficients, originally derived 
by the method of Reed and German (1965), have 
been modified so as to simulate observed tracer 
distributions (agreement with tracer data is discussed 
in the following section). The stream functions 
from which the bulk velocities were calculated were 
based partly on the meridional circulation field 
computed by Cunnold et al. (1975) and partly on 
mean observed meridional winds. 

The model currently contains constituents in the 
oxygen-nitrogen-hydrogen-carbon-chlorine families 
and generally uses the rate coefficients recommended 
in the NASA Chlorofluoromethane Workshop Report 
(Hudson, 1977) or in the work of Turco and Whitten 
(1975) on the effect of chlorofluoromethanes on 
ozone. Exceptions to the rule are discussed later. 

Aerospace Corporation Model 

The Aerospace model (Widhopf, 1975; Widhopf 
et al., 1977) employs spherical geometry, as does the 
Ames model, but it extends from pole to pole in 10° 
intervals; the vertical coordinate extends from the 
surface to an altitude of 50 km with a 1-km grid 
spacing from the surface to 35 km and with a 2.5-km 
grid spacing at higher altitudes. The lateral boundary 
conditions for all constituents are zero flux across the 
polar regions. The species 0(^P), 0(‘D), O3, OH, 
HO2, N, and H are taken to be in photochemical 
equilibrium at the upper boundary; the mass mixing 
ratios of NO2, N2O, CH4, CO, and HNO3 were 
analytically continued to the upper boundary by 
a second-order extrapolation in space and time 
(Widhopf, 1975; Widhopf and Taylor, 1974). At 
the lower boundary the concentrations of O3 , N2 O, 
and CH4 were held fixed (6X10^^ parts/cm 
0.31 ppmv, and 1.35 ppmv, respectively), while 
removal of H2O, HNO3, NO2, NO, and H2O2 from 
the troposphere was simulated by rainout-washout 
(Glatt and Widhopf, 1978; Widhopf and Glatt, 
1979a); H2O and HNO3 were removed at average 
rates defined by Junge (1963) and NO2, NO, and 
H2O2 were assumed to be removed at one-tenth that 
rate. Other species [0(^P), 0(^D), OH, HO2, N, 


and H] were taken to be in photochemical equilib- 
rium at the lower boundary because of their rela- 
tively short lifetimes. Anthropogenic emissions of 
NOx are included in the model, and the lower 
boundary condition for CO is interpreted from 
the measurements of Seiler (1974). 

The chemical rate equations were solved using 
a second-order implicit technique (Widhopf and 
Victoria, 1973). The chemical production-loss term 
Q‘ is approximated by the expansion 



In equation (8) c represents the mass mixing ratio, 
p is the air density, and T is the air temperature. 
Note that equations (7) and (8) are similar to (5) 
and (6), except that (6) does not contain the density 
and temperature terms and is a first-order approxima- 
tion, whereas (7) is a second-order approximation 
in time. Diurnal averaged local photolysis rates are 
computed at every third time step using a technique 
developed by Kramer and Widhopf (1978), and 
diurnal averaging of reaction rate coefficients is 
carried out using a technique similar to that of 
Turco and Whitten (1978). Advective and diffusive 
terms that are important in determining the time- 
dependent distributions of the species are treated 
using a leap-frog and Dufort-Frankel second-order 
accurate finite-difference scheme, respectively. The 
longer lived species (time constant >2 days) are 
treated in an explicit manner. 

The atmospheric meridional bulk velocity and 
eddy transport fields are derived in the manner 
discussed in Widhopf et al. (1977). Briefly stated, 
the mean meridional circulation was obtained from 
the work of Louis (1973) and Louis et al. (1974), 
and the eddy coefficients were based principally on 
the work of Luther (1973), Hunten (1975), and 
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Newell et al. (1966). In testing these coefficients 
against dispersion of inert tracers in the atmosphere, 
they were found not to be totally adequate (Widhopf, 
1975) and were improved by numerical experimenta- 
tion described by Widhopf et al. (1977). 

The important features of the Ames and Aero- 
space models are compared in table 2. 

Reaction Rate Coefficients 

The reaction rate coefficients used in the work 
reported here are essentially those reported in Part I 
with some exceptions; the values are listed in 
tables 3-7. As indicated in those tables, the Aerospace 
model contained no chlorine constituents for this 
particular investigation; although it has since been 
modified to include chlorine, the changes were 
not made in time for this study. Finally, we mention 
the recent measurement (Zahniser and Howard, 
1978) of the rate coefficient for the very important 
reaction (2), the temperature-dependent value of 
which is given as 

A: = (1.4±0.4)X10"^^ 

X exp(-580 ± lOO/T^cm^ sec"^ 

Although this measurement was announced too late 
for most of our assessments, we were able to make 
several calculations with the Ames model using the 
new data for reaction (2). The results are discussed 
in section V. 


IV. AMBIENT ATMOSPHERE 


Before discussing the application of our models to 
impact assessments, it is important to demonstrate 
their capabilities for predicting ambient stratospheric 
characteristics. First, however, we vwU make a few 
remarks about photochemical “time constants.”^ 
If the characteristic time constants are very small, 
a few days or less, the constituent in question is 
nearly in a photochemical steady state (in daytime 
at least); therefore, it is not directly influenced 
by transport, although transport may be indirectly 
important through its influence on the distribution 


3 By “time constants,” we mean the time required for 
the constituent concentration to decrease to of its value 

at the time the source strength is set equal to zero. 


of the parent gas from which the species is formed. 
On the other hand, if the time constants are large, 
abundances and distributions are governed largely by 
transport, and photochemistry may be relatively 
unimportant. An extreme example is carbon 14 — it 
is chemically inert, but does decay radioactively with 
a lifetime in excess of 5,000 yr. Other, less extreme 
examples are the total concentration of odd nitrogen 
(NOy), which is formed within the stratosphere 
and removed, mainly by downward transport, into 
the troposphere, where it can be removed by rainout- 
washout. Other examples are the chlorofluoro- 
methanes (CFM), which are transported from the 
troposphere into the stratosphere where they are 
photolyzed. Because constituents with large time 
constants (of the order of months to years) act as 
tracers, they are especially valuable in studying 
transport in the upper atmosphere; we will now 
consider them in some detail. 


Tracers 

During large thermonuclear explosions, vast quan- 
tities of carbon 14, formed by neutron bombardment 
of nitrogen 14, are injected into the stratosphere 
over a range of altitudes. One can estimate the rates 
of vertical transport by measuring the temporal 
change in the distribution of the (excess) carbon 14 
over a period of several years (Johnston et al., 1976). 
In fact, substantial use has been made of the data 
from the large (57 MT) Soviet high-latitude (75^ N) 
test. However, there may be limitations to the use of 
these data for transport studies: Chang (1975) has 
conjectured that an unsampled reservoir of excess 
carbon 14 was feeding the region at lat. 30° N, and 
Mahlman (1975) showed that lack of a tropospheric 
sink does influence the apparent vertical transport, 
leading to an underestimate of the eddy diffusivity. 
Therefore, one can interpret these simulations as does 
Oliver et al. (1977) as an estimate of the lower bound 
to the rate of vertical transport in the stratosphere. 

Figure 2 shows the predictions of excess carbon 14 
for the two models at lat. 30° N 9 months after 
injection together with the corresponding measured 
number density distributions. The Aerospace model 
is in quite good agreement with observation, but 
the Ames model shows somewhat faster removal of 
the carbon 14 than measurements seem to indicate. 
Here it should be remarked that the Aerospace model 
used the carbon-14 data to refine the eddy transport 
parameterization at higher latitudes. Since the 
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measurements are not extensive enough to provide a 
definitive test of model transport and could be 
interpreted as yielding a lower bound to the rate of 
vertical transport, neither model is inconsistent 
with observation. 

Bauer et al. (1978) have suggested the use of 
zirconium-95 data from recent Chinese atmospheric 
thermonuclear explosions as a tracer of vertical 
transport. However, as they point out, the particulate 
nature of that isotope requires the proper treatment 
of sedimentation velocities. Since neither model is 
able to satisfactorily treat the size distribution 
of such aerosol particles, we do not employ the 
zirconium -95 data in the studies described in this 
report. 

A potentially more definitive set of transport 
tracers (because sources and sinks are known or can 
be determined) is composed of those constituents 
that are transported upward from the troposphere 
into the stratosphere and decomposed there. These 
constituents include N 2 O, CH4, and chlorofluoro- 
methanes. For example, nitrous oxide is a good 
tracer because it has a stratospheric lifetime against 
chemical attack that is very much larger than its 
lifetime against photolysis by sunlight. Hence, using 
knowledge of the photoabsorption cross sections 
and solar radiation intensity to compute the chemical 
loss rates, one can employ such tracers to make 
useful estimates of transport rates in the lower strato- 
sphere. Their usefulness as tracers is complicated 
somewhat by uncertainty as to the location of the 
principal ports of entry into the stratosphere (i.e., are 
they injected primarily in the intertropical con- 
vergence zone (ITCZ) or at other latitudes as well?) 
and lack of adequate data above 30 km. The use of 
methane data as a check of transport parameteriza- 
tion is further complicated by the dominant loss 
mechanism, reaction with OH, because the OH 
abundance is itself dependent on the methane abun- 
dance (see Part I for a complete discussion) and is 
very uncertain. 

Predictions of the vertical distributions of N 2 O are 
compared with measurements in figure 3. Although 
neither model is inconsistent with observations, there 
are insufficient data for a definitive evaluation of 
the transport simulation. The considerable scatter 
in the high-latitude data of Schmeltekopf et al. 
(1977) probably reflects variable conditions of 
transport. Similar considerations hold for the pre- 
dicted and measured methane profiles shown in 
figure 4. This kind of comparison is further com- 
plicated by adoption of a lower boundary condition 


for methane. The tropospheric mixing ratio adopted 
as the lower boundary condition for the Aerospace 
model was 1.35 ppbv, based on early measurements, 
whereas the Ames model employed mixing ratios 
of « 1.8 to 2 ppbv, based on fits to the tropospheric 
and lower stratospheric measurements. More recent 
measurements (Fabian et al. 1979) yield an average 
mixing ratio of approximately 1.65 ppbv. 

To evaluate the simulation of the horizontal as 
well as vertical transport, we turn our attention to 
the predicted distributions of some of the species 
that are formed in whole or in part vrithin the strato- 
sphere: ozone, nitric acid, and water vapor. Ozone 
is formed mainly at altitudes above 30 km due to 
the photodissociation of molecular oxygen in the 
Herzberg continuum (e.g., Whitten and Poppoff, 
1971) and subsequent attachment of the oxygen 
atoms to O 2 . After formation it can, of course, 
be destroyed by the various processes discussed in 
Part I or it can be transported downward and pole- 
ward by atmospheric motions of varying scales. The 
horizontal gradients in the ozone concentration, 
which are associated with the poleward transport, 
are believed to be caused by eddies in the lower 
stratosphere; the eddies, according to Newell (1963), 
cause the mixing surfaces to slope more steeply 
toward the pole than do the surfaces of equal 
potential temperature. That is, parcels of ozone 
moving poleward in the lower stratosphere tend 
also to move downward in the transport region to 
altitudes where their destruction becomes less effi- 
cient. Figure 5 shows some isopleths of ozone con- 
centration for the two models obtained using the 
older value (see Hudson, 1977, for reaction (2)). 
Note that the ozone abundances in the lower strato- 
sphere increase toward the poles and that the altitude 
of peak concentration slopes downward as one goes 
toward the pole. The agreement between the two 
models is quite good except at altitudes below 
~17 km and at high latitudes; the differences below 
17 km may result, in part, from different choices of 
lower boundary conditions for ozone, and those 
at high latitude may reflect different values of the 
transport parameters. 

Figure 6 shows computed and observed ozone 
profiles for three latitudes for the older value of 
the rate coefficient for reaction (2). The ozone 
abundances are overestimated by both models at all 
latitudes. However, figure 7 shows that when the 
newer value for the rate constant (Zahniser and 
Howard, 1978) is used, the agreement between 
observation and theory is greatly improved. The 
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principal reason for the improvement is the substan- 
tially increased ozone loss rate below 25 km (see 
Whitten et al,, 1978; Turco et al., 1978). Figure 8 
shows the corresponding computed and observed 
column densities for the summer and faU seasons; 
obviously, the agreement with observation is much 
improved. Similar results are obtained in other 
studies by Widhopf and Glatt (1979b). Here it 
should be remarked that neither model contained 
chlorine when the computations with the older 
value of the rate coefficient for reaction (2) was 
used. However, the Ames model was altered to 
simulate the effect of 2.5 ppbv of Clx (the sum of 
HCl, Cl, CIO, and CIONO 2 ). The absence in the 
simulation of observed abundances of chlorine com- 
pounds could not account for the overestimate of 
stratospheric ozone abundance. 

A second tracer that is formed within the strato- 
sphere is “odd-nitrogen,” composed of NO, NO2, 
NO3, N2O5, HNO3, HNO2, and CIONO2 and some- 
times called “NOy” (the sum of NO and NO2 con- 
centrations is usually called “NOx”)* Although 
odd-nitrogen is not a directly measured quantity, 
its principal components, NO, NO2, and HNO3, 
have been measured individually. Total NOy is 
important for an intercomparison of the two models 
because (1) the comparison of its predicted distri- 
bution shows the effects of the two different trans- 
port parameterizations, and (2) the predicted amount 
and distribution of NOy in the atmosphere is related 
to the predicted change in ozone. Figure 9 shows 
profiles of computed NOy at low, mid-, and high 
latitudes and some (mostly midlatitude) data repre- 
sentative of a compilation by Ackerman (1975) of 
numerous measurements of NO, NO2, and HNO3. 
At altitudes between 17 and 40 km, the predictions 
obtained from the two models agree with each other 
quite well. However, at low altitudes the Aerospace 
model predicts substantially less NOy than does the 
Ames model, except in the lower 4-5 km where 
anthropogenic NOx injection is simulated in the 
Aerospace model. At high altitudes, the Aerospace 
model predicts more NOy than does the Ames model 
because the former assumes a small source of meso- 
spheric NOy that is absent from the latter. Such a 
source may be characteristic of auroral latitudes 
where electron bombardment dissociation of N 2 is 
expected to yield a net production of NO. At lower 
latitudes (i.e., equatorward of '--60°) there probably 
exists a net NO sink through the reaction 

N + NO ^ N 2 + O (9) 


This difference in upper boundary conditions effects 
the distributions appreciably only above ~35 km 
(Duewer et al., 1977). The differences at altitudes 
between 14 and 17 km are due, in part, to the simula- 
tions of washout-rainout in the Aerospace model; 
washout-rainout is not simulated in the Ames model. 
We believe these differences to be useful in exploring 
the model dependence of ozone destruction discussed 
in section V because they provide limiting cases for 
natural sources and sinks. 

The principal NOy species in the lower strato- 
sphere is nitric acid; since many observations of this 
constituent have been made (Lazrus and Ganrud, 
1974; Murcray et al., 1975), it is useful to compare 
predictions with the observations. However, column 
densities show much less variability due to small- 
scale variations in transport than do local concentra- 
tions. Therefore, the column densities that cover a 
range of altitudes and show much less variability 
are most useful for comparisons with model predic- 
tions. Such comparisons are shown in figure 10. 
Obviously the predicted and measured column 
densities for the two models greatly exceed the 
measurements, although the predicted latitudinal 
gradients agree rather well with the measured 
gradients. 

In order to understand the discrepancy, it is 
useful to examine another related parameter, the 
HNO3/NO2 ratio, and to compare the calculated and 
observed values. The usefulness of the comparison 
is apparent from considerations of the NOy budget. 
If the nitric acid abundance is high compared with 
observation, one would expect the HNO3/NO2 
ratio to be large. Estimates of the ratio as a function 
of altitude obtained by several investigators are 
shown in figure 1 1 together with predictions obtained 
with the Ames model for lat. 45° N and with the 
Aerospace model at lat. 40° N. Although the model 
predictions strongly disagree with the data of Evans 
et al. (1976), they do not disagree nearly as strongly 
with the other data. In any case, the high predicted 
HNO3 levels shown in figure 10 are roughly consis- 
tent with the apparently high predicted HNO3/NO2 
ratios shown in figure 11. The Aerospace values are 
higher than the Ames ratios, which is consistent 
with the larger HNO3 column densities obtained 
in the Aerospace model. Schiff et al. (1978) and 
Evans et al. (1978) have suggested that the discrep- 
ancy between measurement and theory can be 
resolved by assigning a pressure dependence to 
reaction (1) such that the rate coefficient is substan- 
tially reduced; however, at present it is not clear 
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that such a pressure dependence exists.^ The relation 
of the HNO 3 /NO 2 ratio to predicted effects on 
stratospheric ozone of SST-emitted NOx will be 
discussed in section V. 

A third species which has many of the attributes 
of a tracer is stratospheric water vapor. Unfortu- 
nately, the sources of water vapor are poorly under- 
stood. Many years ago Dobson et al. (1946) proposed 
that water vapor enters the stratosphere almost 
entirely through the tropical tropopause, which forms 
a cold trap; the cold trap then determines the 
water-vapor mixing ratio. Since then, it has become 
generally accepted that methane oxidation yields 
an additional source, largely in the middle strato- 
sphere. Harries (1976) has summarized the various 
measurements made up to 1976; the range of mass 
mixing ratio values, obtained by averaging the mea- 
surements, is shown in figure 12. Comparison with 
the predicted mixing ratios, also shown in figure 1 2 , 
indicates that the agreement is reasonable; however, 
the wide spread of the measured values should be 
noted. Note that the Aerospace model calculates the 
distribution of water vapor everywhere and that the 
Ames model prescribes the tropospheric distribution 
using standard atmosphere data. 


Reactive Species 

There are several highly reactive species that are 
important to our studies: OH, HO 2 , NO, N 02 ,and 
CH 3 O 2 . We will begin our discussion with the odd- 
hydrogen radicals OH and HO 2 . 

The hydroxyl radical is formed in the stratosphere 
by the reaction of excited oxygen atoms with water 
vapor: 

0(^D) + H 2 0 20H (10) 

and by methane oxidation. The latter process, which 
is initiated by the OH destruction reaction 

OH + CH 4 ^ H 2 O + CH 3 (11) 

is believed to yield a net production of nearly two 
OH molecules for each methane decomposition 
reaction (see Part I for a complete discussion). 
In the upper stratosphere the methane oxidation 


^ If this reaction of OH with HNO3 should prove to be 
substantially faster than indicated in table 2 , the HNO3/NO2 
ratio could be significantly lowered. 


sequence is a small contributor to OH formation, 
but in the lower stratosphere, as our studies have 
shown, it is nearly as important as the water-vapor 
source of OH. The hydrogen radicals are lost pri- 
marily by reaction with each other, 

OH + HO2 H2O + O2 ( 12 ) 

We have plotted the predictions of OH concen- 
tration in figure 13, which also shows Anderson’s 
observed values obtained at a solar zenith angle of 
^^ 86 ° (at Palestine, Texas, lat. 32° N, in July and 
January). Because of the large solar zenith angle, 
the observed values should be multiplied by a factor 
of about 1.5 to 2 when compared with the mean 
daytime predicted concentrations shown in fig- 
ure 13(a). Hence, they appear to be somewhat 
larger than our predictions would indicate for large 
solar zenith angles, at least for heights above 30 km. 
Furthermore, we note that Anderson’s January 1976 
values were obtained a little later in the year than 
were the corresponding predicted values (mid- 
November) and at a slightly lower latitude; in mid- 
February at lat. 35° N, the predicted concentrations 
of OH are about 20% lower than those shown for 
lat. 40° N in mid-November. Thus, the agreement 
between measured and predicted hydroxyl con- 
centrations can be regarded only as fair for altitudes 
above 30 km, although the difference between 
observation and calculation is much less than it was 
prior to the measurement of the rate coefficient for 
reaction ( 1 ) by Howard and Evenson (1977). 

So far, no stratospheric measurements of OH 
abundance have been made below 30 km. At altitudes 
below 20 km, there are substantial differences 
between the predictions of the two models relative 
to OH, which we attribute largely to the differences 
in predicted ozone concentrations and rainout- 
washout of HOx and HNO3. Smaller ozone abun- 
dances, as obtained with the Aerospace model, lead 
to less predicted 0(^D) and thus a small rate of 
formation of OH. The absence of rainout-washout 
in the Ames model results in higher HOx, NOx, 
and HNO3 abundances. Comparison with HNO3 
measurements and NOx estimates due to Fishman 
and Crutzen, obtained by balancing the CO budget, 
show good agreement with Aerospace model predic- 
tions in the troposphere (Widhopf and Glatt, 1979a). 
Smaller OH concentrations also mean that the rate 
of OH conversion to HO2 via the reaction 

OH + O3 HO2+O2 ( 13 ) 
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is slowed, thus decreasing the effect on OH of the 
smaller ozone abundance. The larger OH concen- 
trations predicted by the Ames model at altitudes 
below 20 km have an important bearing on the 
differences in the SST effects on ozone predicted 
by the two models. This aspect is discussed in 
section V, 

Figure 14 compares predictions of HO2 for the 
two models. Although no measurements of HO2 
have yet been made, we show these results as a 
further useful comparison of the HOx chemistry 
contained in the models. We attribute the substantial 
difference at high latitudes in autumn to differences 
in computed mean solar zenith angles near lat. 70° N. 
The radical CH3O2 , which has not yet been observed 
in the stratosphere, can also be classed as an HOx 
radical because of the similar chemical behavior of 
H and CH3. The predicted profiles of CH3O2 are 
substantially different above 30 km, as shown in 
figure 15; however, it is not important to the 
chemical cycle above that altitude. The reason for 
this divergence is the dependence of CH3O2 produc- 
tion on the distribution of its methane source 
(cf. fig. 4). Because it is less abundant and probably 
somewhat less reactive, the radical CH3O2 is 
probably of much less direct importance to ozone 
chemistry than is HO2. Its main role is as an inter- 
mediary in the methane oxidation sequence. 

As has been discussed many times (see Part I for 
references) the odd-nitrogen species NO and NO2 
strongly influence ambient ozone as well as ozone 
loss due to supersonic transport emissions. Figure 16 
compares the distributions predicted by the two 
models with each other and with the results of a 
simultaneous measurement of NO and NO2 by 
Ackerman et al. (1975). It is seen that the predictions 
of both models are slightly smaller than the measured 
values, but the disagreement is not serious (see the 
HNO3/NO2 data in fig. 11). The principal reasons 
for the differences between the two sets of predic- 
tions are the effects of the tropospheric conditions 
(discussed earlier with respect to NOy) and slightly 
different partitioning of NOx between NO and NO2- 
Partitioning is determined approximately by the 
relation 

(14) 

[NO 2 ] A:[03] ^ ^ 


where J is the photolysis rate of NO2 , [O3 ] is the 
local ozone number density, and k is the rate coeffi- 
cient of the reaction 

NO + O3 NO2+O2 (15) 

If partitioning favors NO2 , more of the NOx will be 
converted to nitric acid and total NOx will decrease 
slightly. The partitioning of NOx predicted by the 
two models is very close at altitudes above 20 km; 
at lower altitudes, the Aerospace model predicts 
less NO2 and more NO than the Ames model because 
of the smaller ozone abundance computed by the 
former. 


Conclusions 

We can summarize the principal points made in the 
preceding discussion of the ambient (unperturbed) 
atmosphere as follows: 

1. Insufficient data are available on the distribu- 
tions of excess carbon 14 from atmospheric nuclear 
tests and on the ambient distributions of nitrous 
oxide and methane to permit the transport param- 
eters of one model to be chosen over those of the 
other. Although the ozone and NOy data appear 
to be in reasonable agreement with observations, we 
must nevertheless conclude that insufficient tracer 
data are available for a satisfactory definition of the 
transport parameters. 

2. The predicted ozone abundance computed by 
both models is in reasonable agreement with observa- 
tion when the new rate coefficients for reactions (1) 
and (2) are used. 

3. The NOy computed by the two models agrees 
quite well between 18 and 40 km altitude. The 
differences below 1 8 km and above 40 km are due to 
different boundary conditions and rainout-washout. 

4. The predicted OH abundance at altitudes 
above 30 km appears to be low in light of Anderson’s 
(1976) measured values. On the other hand, it may 
be too high at lower altitudes, which is suggested 
by the apparently excessive computed nitric acid 
concentrations (fig. 10) and HNO3/NO2 ratios 
(fig. 11). The latter characteristics are important 
because they are closely related to the magnitudes 
of the predicted ozone changes induced by NOx 
emissions from SST’s. 
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V. RESULTS AND DISCUSSION 


In this section we present and discuss the results 
obtained for operation of types A and B aircraft 
(20-km and 17.5-km nominal flight altitudes, 
respectively), using the emissions model shown in 
figure 1 ; one scenario for the type A aircraft includes 
emissions by subsonic aircraft as well. These results 
take the form of plots of relative ozone column 
change as a function of latitude and the vertical 
distributions of the absolute ozone and NOy changes 
at various latitudes. We then compare the results from 
the Ames and Aerospace models with each other 
and with the results of previous multidimensional 
model calculations (e.g., Alyea et al., 1975; Borucki 
et al., 1976; Widhopf et al., 1977; and Hidalgo and 
Crutzen, 1977). Finally, we briefly discuss the 
uncertainties associated with the calculations. 


Flight Scenarios 

To arrive at reasonable estimates of the possible 
effects on stratospheric ozone of supersonic aircraft 
operations as described in section II, we constructed 
the following scenarios: 

1. Five years of operation of type A aircraft 
(nominal cruising altitude of 20 km), using the NOx 
injection rate given in figure 1 and the corresponding 
water-vapor injection rate (~200 g H 2 0/g NO 2 ). 
In the Ames model, the simulation of NOx 
water-vapor emission was performed at an altitude 
of 20 km because of the 2.5-km vertical grid spacing. 
For the Aerospace model, injection at 19, 20, and 
21 km was simulated, as shown in figure 1. No 
chlorine constituents were included in the models 
and the older value of the rate coefficient for reac- 
tion (2) was used. We note here that when the new 
value for that rate coefficient was announced by 
Zahniser and Howard (1978), the 5-yr runs had 
already been completed. 

2. Five years of operation of type B aircraft 
(nominal cruising altitude of 17.5 km). The altitudes 
at which injection was simulated were lowered: 
to 17.5 km for the Ames model, and to 16, 17, and 
18 km for the Aerospace model. Otherwise, the NOx 
emission rates shown in figure 1, along with the 
corresponding water-vapor emission rates, were 
used. The same chemical model used to perform 
the case of injection at an altitude of 20 km was 


also used here because these runs were also completed 
prior to the announcement of the new rate for 
reaction (2). 

3. Five years of operation of type A aircraft 

together with a nominal fleet of subsonic aircraft, 
using the emission rates shown in figure 1 and table 1. 
To simulate the removal of NOx troposphere, 

the increase in HNO3 was set equal to zero at an 
altitude of 5 km in the Ames model. A comparison 
of our results for the increase of NOy (ANOy) with 
results from a variant of the Ames model which 
contained “rainout” boundary conditions showed 
that causing ANOy to vanish at an altitude of 5 km 
gave a good approximation to the “rainout” condi- 
tions. The Aerospace model simulated rainout- 
washout directly in the same manner discussed for 
the calculation of the natural atmosphere. 

In order to then simulate steady-state conditions 
in an atmosphere containing chlorine constituents, 
we employed two additional scenarios, using the 
Ames model only : 

4. Fifteen years of operation of type A aircraft 
using the NOx injection rate given in figure 1; to 
minimize computer usage, the NOy injection, without 
any chemistry, was simulated for 13 yr; this simplifi- 
cation is justified because stratospheric NOy is a 
nearly-conserved quantity. The full model, which 
contained 2.5 ppbv of ozone-active chlorine (at an 
altitude of 60 km), was then run until the seasonal 
changes in ozone abundance became repetitive 
(~2 yr); the last two seasons (summer and autumn) 
were then very close to a steady-state cyclic variation. 
The new value of the rate coefficient of reaction (2) 
was used. 

5. Same as (4), but with type B aircraft. 


Changes in NOy 

Figures 17 through 21 show the changes in NOy 
(ANOy) that are predicted by the models for 
scenarios (1) to (5), respectively, at latitudes of 0°, 
40°, and 70° N. We present results for only the 
Northern Hemisphere because nearly all of the 
simulated aircraft emissions occur there; the eventual 
buildup of ANOy in the Northern Hemisphere would 
be nearly twice as large as in the Southern Hemi- 
sphere. Plots of ANOy are useful because the distri- 
bution is influenced only by transport and the 
mesospheric and tropospheric sinks; stratospheric 
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photochemistry does not affect it. The predicted 
spatial distributions of the increase in water vapor 
(AH 2 O) due to SST emissions was found to be very 
similar to those of ANOy (except for a multiplicative 
factor of about 550); hence, we do not show them 
here. 

Comparison of the vertical profiles obtained by 
the Ames and Aerospace models for scenario (1), 
shown in figure 17, and scenario (2), presented in 
figure 18, suggest that the Aerospace model predicts 
that ANOy is concentrated in a narrow altitude 
range. This is a result of at least three effects. First, 
the vertical resolution is quite different (1 km vs 
2.5 km), which will tend to maintain a sharper 
profile in the Aerospace model. Second, the Aero- 
space model includes rainout-washout, which tends to 
decrease the concentration in the troposphere faster. 
Finally, the faster vertical diffusion in the Ames 
model accounts for some of the difference. Such 
relative behavior of the predicted ANOy is consistent 
with the predicted temporal behavior of excess 
carbon 14 shown in figure 2. In addition to the 
differences in vertical diffusion, differences in resolu- 
tion and transport evidently cause the maxima in 
the predicted ANOy profiles to occur at an altitude 
that is about 2 km lower at lat. 70° N in the 
Aerospace model than in the Ajnes model. As we 
noted previously, neither model is inconsistent with 
available tracer observations; hence, the results 
from each model should be regarded as reasonable 
simulations of the distribution of the ANOy. 

The distribution of the ANOy is rather strongly 
dependent on the flight altitude. For the case of 
NOy injection at 20 km, more of the simulated 
NOy is retained in the stratosphere near the injection 
altitude by the Aerospace model than by the Ames 
model (fig, 17). However, the Ames model predicts 
faster upward transport to higher altitudes for 
injection at 17.5 km. At high latitudes the Aerospace 
model tends to retain more of the ANOy at lower 
altitudes than does the Ames model, and this has a 
significant bearing on the ozone change resulting 
from stratospheric injection. This is probably affected 
by washout-rainout, which would decrease ANOy 
proportionately faster for the 1 7.5-km case than for 
the 20-km case, and by the effect of faster vertical 
diffusion in the Ames model. In the 17.5-km- 
injection-height case more total ANOy is retained by 
the Ames model than by the Aerospace model, 
which is the opposite of the carbon- 14 result. This 
indicates that rainout-washout accounts for the 
differences, because such processes are not important 


to carbon-14 distributions. As we shall see later 
(also see Part I), the increase of NOy at relatively 
low altitudes (i.e., below -^20 km) leads to ozone 
increases, and enhanced NOy at higher altitudes 
leads to ozone decreases. One should also note in 
both models the downward trend of the ANOy as 
it is moved poleward. Such an effect is, of course, 
expected because it is closely related to two high- 
latitude phenomena: (1) the increase in ambient 
ozone and (2) the lower altitudes at which peak 
ozone concentrations occur. That is, some of its 
sources and sinks are similar to those of ozone in 
that NOy is formed in the upper stratosphere (via 
reaction of 0(^D) with N 2 O) and lost by transport 
into the troposphere. 

In scenario (3) the Aerospace model, as we might 
expect from our foregoing comments concerning 
vertical diffusion, resolution, and rainout-washout, 
shows retention of the vertical distribution of 
injected NOx (see table 1) and less vertical diffusion 
than the Ames model. In this the former predicts 
a deeper minimum between the stratospheric and 
tropospheric contributions at mid-latitudes than 
does the Ames model. 

We cannot, of course, compare results from both 
models for scenarios (4) and (5) (figs. 20 and 21); 
however, it is worth noting that the vertical distri- 
bution of the injected NOy is more dispersed in the 
steady-state case than after 5 yr of operations. The 
spreading is asymmetric; that is, the ANOy increase 
is greater above the peak than below it. From our 
comments above, such behavior has an important 
bearing on predicted ozone changes since the more 
NOy that is transported to higher altitudes, the 
greater the likelihood that ozone will decrease. 

Changes in Ozone 

Predicted changes in the vertical distribution of 
ozone (AO 3 ) as a consequence of stratospheric 
emissions of NOy and water vapor are presented in 
figures 22 through 27. The most striking features 
of the profiles of AO 3 are the changes from ozone 
gains to ozone losses at altitudes between 18 and 
22 km. In general, calculations using the Ames model 
showed smaller ozone gains at low altitudes and larger 
ozone losses at high altitudes than did the Aerospace 
model. Such behavior of the stratospheric ozone is 
in part a consequence of the predicted vertical dis- 
tributions of ANOy shown in figures 17 through 21 
and in part a consequence of the vertical distributions 
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of OH and HO2. Note that at altitudes between 
10 and 20 km the Ames model contains more OH 
than the Aerospace model, and between 20 and 
30 km, less HO2 than the Aerospace model (except 
at high latitudes in autumn). The overall effect of 
these differences is a prediction by the Aerospace 
model of more ozone gain between 10 and 30 km. 
Figure 24 includes a simulated fleet of subsonic 
aircraft (see table I) along with a nominal fleet of 
type A SST’s. In this case the ozone gains, which 
occur in the troposphere, greatly exceed the losses, 
which occur above ~20 km, and thus yield a net 
increase in ozone. For scenarios (1) and (2) the 
changes in total ozone column are very close to zero 
for the Aerospace model, but are negative and of the 
order of tenths of a percent for the Ames model. 
Discussion of the Ames predictions of ozone column 
change is deferred until we discuss the steady -state 
results. 

The catalytic chain of reactions responsible for 
ozone destruction by nitrogen oxides has been 
known for a long time (Crutzen, 1970; Johnston, 
1971): 


NO + 03 

^ NO2 + O 2 

(16) 

NO2 +0 

^ NO + O, 

(17) 

Net: O + O3 

2O2 


and is not discussed 
writing the catalytic 

further here (see Part I). In 
sequence composed of reac- 


tions (16) and (17), one must keep in mind that the 
quantity of NO 2 participating in reaction (17) is 
exactly that produced in reaction (16). Other sources 
of NO2 must be excluded when evaluating the net 
rate of the cycle. Similar arguments hold for the 
other cycles discussed below. Because of reac- 
tions (16) and (17), increases of NOy at higher 
stratospheric altitudes lead to decreases in ozone. On 
the other hand, increases of NOy at lower altitudes 
lead to increases in ozone because of reactions that 
have only recently been recognized; they are reac- 
tions that catalytically form ozone (see Part I; also 
Widhopf et al., 1977; Hidalgo and Crutzen, 1977; 
Crutzen and Howard, 1978; Turco et al., 1978): 


NO + HO 2 

-> 

NOj + OH 

(1) 

^ NO 2 + hp 


NO + 0 

(18) 

0 + O 2 + M 


O 3 +M 

(19) 

Net: HO 2 + O 2 


OH + O 3 



and to a lesser extent 
NO + CH3O2 


NO2 +CH3O 

(20) 

NO2 +hn 

- 

NO + O 

( 18 ) 

0 + O2 + M 


O3+M 

( 19 ) 

Net: CH3O2 + O2 


CH3O + O3 



in which CH3O2 is an intermediate product in the 
methane oxidation chain. Although reaction (20) 
provides a direct link between methane oxidation 
and ozone generation, a more important, but 
indirect, influence results from the copious forma- 
tion of HO2 during the methane oxidation process. 
In fact, methane decomposition may be the largest 
source of HO2 in the lower stratosphere. In addition 
to ozone formation throu^ the sequences listed 
above, NO introduced into the lower stratosphere 
converts HO2 normally present into OH. Such 
conversion has two effects: (1) it shifts the reactive 
NO component of NOy into the relatively inert nitric 
acid and (2) it amplifies the oxidation rate of 
methane and CO, thereby generating HO2 and ozone 
more rapidly. 

As discussed elsewhere (Part I; Crutzen and 
Howard, 1978; Turco et al., 1978) the earlier failure 
to recognize the significance of reactions (1) and (2) 
was not recognized due to the low values ascribed 
to the rate coefficients (see Hudson, 1977, for a 
discussion relative to reaction (1), and Whitten 
et al., 1978, for a discussion relative to reaction (2)). 
The larger, and now generally accepted, values 
(see table 2) were first reported by Howard and 
co-workers (Howard and Evenson, 1977; Zahniser 
and Howard, 1978). 

In computing the ozone changes shown in fig- 
ures 25 and 26 (i.e., steady -state conditions), an 
ozone-active chlorine mixing ratio of 2.5 ppbv at an 
altitude of 60 km was included in the Ames model 
because such a chlorine abundance is probably a 
reasonable extrapolation to the years 1990-2000. 
Chlorine is important in the present context because 
it is closely coupled to NOx through the reactions 

CIO + NO2 + M ^ CIONO2 + M (21) 

CIO + NO ^ CI + NO2 (22) 

Thus, if stratospheric NOx is increased, the rates of 
reactions (21) and (22) will be increased. Since reac- 
tion (21) tends to shift Clx away from the active 
forms, Cl and CIO, increasing its rate will decrease 
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ozone loss due to chlorine catalysis. Reaction (22) 
competes with the rate-controlling odd-oxygen loss 
process 

ClO + 0 ^ CI + O2 (23) 

Hence, increasing the rate of reaction (22) will 
decrease the effectiveness of reaction (23). As we 
discussed in Part I, the emission of NOx at SST 
flight altitudes can therefore be expected to lead to 
greater ozone increases (or smaller ozone losses) 
when substantial abundances of chlorine are present. 

As noted before, our studies also included the 
simulation of water-vapor emission by aircraft flying 
in the stratosphere. Added water vapor is important 
because, through its oxidation by 0(^D), it leads 
to enhanced HO2 and OH abundances. Although 
increased HO2 does have a tendency to increase 
ozone via reaction (1), such an effect is overwhelmed 


by direct destruction by the HOx^ 

OH + O3 ^ HO2+O2 ( 14 ) 

HO2+Q ^ OH + O2 ( 24 ) 

Net: O + O3 2O2 

HO2 +03^ OH + 2O2 ( 2 ) 

OH + O3 ^ HO2 + O 2 ( 14 ) 

Net: O3 + O3 3O2 


Thus, a large value for the rate coefficient for reac- 
tion (2) can be expected to lead to non-negligible 
ozone decreases (Turco et al., 1978). Because OH 
attacks the dominant stratospheric chlorine reservoir, 
HCl, to release the ozone-active Cl, 

OH + HCl H2O + CI (25) 

an increased stratospheric OH abundance can be 
expected to lead to greater ozone destruction by 
chlorine catalysis. In fact, we showed in Part I, that 
in a globally averaged sense and for the aircraft- 
emission indices considered, HOx - Clx coupling is 
expected to result in an ozone loss that nearly 
balances the ozone gain resulting from NOx - Clx 
coupling. 

It is important to point out here that our 
model calculations did not take account of the 
thermal feedback effect: an ozone decrease will lead 
to lower upper stratospheric temperatures; because 

14 


of substantial positive activation energies of the 
reactions that govern the catalytic destruction of 
ozone, such temperature changes will tend to limit 
ozone loss. For the case of water-vapor injection, 
increased stratospheric water-vapor loading causes 
the cooling rates to increase because of increased 
emission of thermal infrared radiation by H2O. 
Stratospheric temperatures are thus lowered, 
causing the catalytic reactions that destroy ozone 
to slow. Luther and Duewer (1977) carried out 
such computations for a case in which they assumed 
the water-vapor mixing ratio increase to be constant 
with altitude, and found that thermal feedback^nearly 
canceled the destruction that would have occurred if 
such feedback were absent. Of course, as figures 17 
through 21 suggest, the AH2O can be expected to 
be much larger below 30 km than above. Since 
thermal feedback is probably significant only above 
30 km, we expect it to be of less importance than 
was proposed by Luther and Duewer. However, 
Glatt and Widhopf (1979a) studied the temperature 
coupling using a two-dimensional steady-state model 
for a similar combined supersonic and subsonic fleet 
and found heating below the injection region, 
whereas cooling occurred in and above the injection 
region. Thus, temperature feedback may be impor- 
tant in the lower stratosphere. 

To investigate the limiting situation in which 
ozone destruction due to water-vapor emission is 
absent, we simulated a variant of scenario (5) in 
which we omitted the AH2O from the calculations. 
The results are shown in figure 27. As expected, 
the ozone gains below ~25 km are much greater 
than those shown in figure 25 (in which water- 
vapor emission is included), and the ozone losses 
above 25 km are less than indicated in figure 25 
(cf. fig. 13(a) of Part I). 

Figures 28 through 32 show the computed relative 
ozone column changes using the Ames model 
(scenarios (4) and (5)); we assumed steady-state 
conditions in an atmosphere containing a mixing 
ratio of 2.5 ppbv of ozone-active chlorine at an 
altitude of 60 km. For the scenario (4), shown in 
figure 28, we assumed operation of type A aircraft 
(20-km nominal cruising altitude) and the old 
value of the rate coefficient for reaction (2). 
Small ozone decreases due to NOx and H2 O emission 
are predicted at all latitudes and seasons, but their 
magnitudes are somewhat larger in the Northern 
than in the Southern Hemisphere. We can qualita- 
tively account for the results shown in figure 28 
by referring to figure 20, which shows that in summer 
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and autumn ANOy in the catalytically destructive 
region above 20 km is larger at high latitudes 
(Northern Hemisphere) than at lower latitudes. 
Because catalytic destruction dominates ozone 
production under the conditions of figure 28, we 
find that the latitudinal behavior of the change in 
ozone column shown is quite reasonable. To clarify 
our simulation of the interaction of added NOx 
with ozone when the new large value of the rate 
coefficient for reaction (2) is used, the results for 
scenario (4) are presented in figure 29. The results 
shown in figure 29 are for the addition of NOx 
only by a fleet of supersonic aircraft operating at 
a nominal cruise altitude of 20 km; there is a small 
increase in the ozone column at all latitudes. Further- 
more, the ozone increase is nearly independent of 
season and is also nearly symmetric about the 
equator. Figure 30 shows the expected increase as a 
function of latitude when water-vapor emission only 
is included in the computations. Ozone depletion 
caused by water-vapor emission is a minimum near 
the equator (~0.7%) rising to larger values at high 
latitudes (~1.4% above lat. 60° N, '--1.1% below 
lat. 60° S. The difference between the corresponding 
abscissas of figures 29 and 30 gives (approximately) 
the predicted ozone change due to NOx water- 
vapor emissions, which is presented in figure 3 1 . The 
corresponding ozone change for scenario (5) is 
presented in figure 32. The slight hemispherical 
asymmetry in ozone loss due to water-vapor emission 
is largely responsible for the rather strong asymmetry 
shown in figure 31. Both the structure in the 
percentage ozone loss and the latitude-dependence 
portrayed in figures 31 and 32 are quite uncertain; 
the total column loss, which is very close to zero, 
is the remnant of the near-cancellation of ozone loss 
at high altitude and gain at low altitude. 

Comparison with Results of Earlier Studies 

Many one-dimensional model studies using the 
older values of the rate coefficients for reactions (1) 
have been carried out (see Grobecker et al., 1974, 
and Part I for a review of the results). However, to 
our knowledge, the only one-dimensional model 
calculations using the new, larger values of the 
rate coefficients for both reactions (1) and (2) are 
those reported by Turco et al. (1978). It is instructive 
to compare the globally averaged ozone change for 
scenario (4) vrith and without water- vapor emission 


(figs. 30, 31) to the corresponding one-dimensional 
model predictions (fig. 2 of Turco et al., 1978). The 
globally averaged equivalent NOy injection rate is 
0.58X10® kg NO 2 per year; the predicted ozone 
increase for the Chang-Dickinson eddy-diffusivity 
profile and no water-vapor emission is —1.5%, which 
is very close to the value —1.25% given in figure 30. 
For the case that includes the effects of water-vapor 
emission, Turco et al. gave results only for the Wofsy- 
type eddy diffusivity, —0.2% of ozone gain which is 
to be compared with the 0.21% given in figure 28. 

All of the earlier multidimensional studies (Alyea 
et al., 1975; Borucki et al., 1976; Widhopf et al., 
1977; and Hidalgo and Crutzen, 1977) used the old 
(pre-1977) values of the rate coefficients of reac- 
tions (1) and (2); they also neglected the effects of 
water-vapor emissions, which was justifiable in the 
light of the pre-1977 knowledge of stratospheric 
chemistry. As a result, simulations of supersonic 
aircraft operations led to predictions of ozone loss, 
globally averaged relative reductions being roughly 
consistent with predictions from one-dimensional 
models. However, the last two papers reported results 
obtained for subsonic aircraft flying in the tropo- 
sphere and in the low stratosphere (below —14 km). 
The low altitudes of injection of NOx ^^rom subsonic 
aircraft generally was shown to yield predictions of 
ozone increase because most of the NOx was retained 
in the lower atmosphere where it generated ozone via 
the methane oxidation sequence and reaction (20); 
because little of the NOx emitted at low altitude 
rises to the region of catalytic destruction of ozone, 
the net result is an ozone increase. However, Widhopf 
and Glatt (1979a) did include water-vapor modeling 
and found that with the use of the new rate coeffi- 
cient for reaction (1) a predicted increase in both 
stratospheric and tropospheric ozone resulted. In 
all of the studies of NOx emissions from supersonic 
aircraft in which it was assumed that most of the 
emission occurred in northern midlatitudes, much 
larger effects were predicted in the Northern than 
in the Southern Hemisphere. Furthermore, the 
predicted relative ozone decrease increased with 
increasing latitudes (e.g., see fig. 14 of Borucki et al., 
1976, for a comparison of three model results). 
This latitudinal dependence was also obtained in the 
current studies in which the Ames model was used 
when the old value of the rate coefficient for reac- 
tion (2) was used and a nominal cruise altitude 
(type A aircraft) was assumed (see fig. 28), although 
the latitudinal gradient is not nearly so pronounced 
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as when the old value for the rate coefficient for 
reaction (1) was used. The reason for that charac- 
teristic of figure 28 is that catalytic destruction of 
ozone outweighs low-altitude ozone generation. As 
shown in figures 29 to 32, the latitudinal dependence 
of the ozone decrease is quite different when the new 
values of the rate coefficients for reactions (1) and 
(2) are used. The differences occur because the model 
then predicts that more ozone is generated at low 
altitudes by the NOy addition and destroyed at 
high altitudes by the water-vapor addition. Ozone at 
lower altitudes is largely controlled by transport; 
hence, the greater the percentage of total ozone 
at lower altitudes, the greater the influence of 
meridional transport. Including the effect of washout- 
rainout on NOx may reduce the effects seen in the 
Southern Hemisphere, since the lower altitude NOy 
would then be taken out of the troposphere faster 
than it could be transported to the Southern 
Hemisphere. 

Sensitivities and Uncertainties 

Because of the large computer times required 
for two-dimensional model simulations, it is not 
practical to use them for sensitivity studies. However, 
a number of sensitivity tests using a one-dimensional 
model were reported in Part I (table 7); our estimates 
of sensitivity are expected to apply to our two- 
dimensional models (in a globally averaged sense) 
as well. We note here the sensitivity study carried 
out with the rate coefficient for reaction (2) (varying 
its value by a factor of 3); the actual increase of the 
rate coefflcient value by a factor of '^4.5 as a conse- 
quence of its measurement by Zahniser and Howard 
(1978) leads to changes in predicted ozone perturba- 
tion that are qualitatively consistent with the ozone 
change differences between figures 28 and 3 1 . 

Table 8 of Part I summarizes most of the sources 
of uncertainty in our SST ozone perturbation 
estimates. However, it is necessary to qualify our 
discussion of transport uncertainty because we are 
dealing here with horizontal as well as vertical trans- 
port. In our earlier discussion of the distributions of 
trace species and their prediction, we showed that 
although the transport parameters of the Ames and 
Aerospace models are not inconsistent with observa- 
tion, neither are they firmly supported by observa- 
tion. In fact, the substantial differences between the 
transport models, indicated in section IV and in 


figures 17 through 21, show that they cannot both 
be accurate simulations of nature. These differences 
are mainly responsible for the differences in pre- 
dicted ozone perturbations that were obtained from 
the two models: After 5 yr of simulated type A SST 
operations (and using the old value of the rate coeffi- 
cient for reaction (2) the Ames model predicted 
small ozone decreases, and the Aerospace model 
predicted very small ozone increases. In assessing 
the significance of the uncertainty associated with 
transport, one should bear in mind the small 
magnitudes of ozone change that are involved and 
that they result from the near-cancellation effect 
of large values of catalytic destruction of ozone at 
high altitudes pitted against large values of ozone 
generation at low altitudes. Thus, small changes in 
the transport parameters can have large relative 
effects on predicted ozone gain or loss. 

VI. CONCLUDING REMARKS 


From the foregoing discussion we see that not 
only have supersonic transport assessments changed 
since the Climatic Impact Assessment Program, but 
they have even undergone modifications since the 
publication of Part I of the present work (Poppoff 
et al., 1978). Although figure 16 in Part I pertains 
to the history of one-dimensional model predictions, 
it is appropriate to update the figure by remarking 
that the estimate of 0.2 to 0.3% ozone gain opposite 
the “no H 2 O” should be increased. With the new, 
larger rate coefficient for reaction (2) the estimates 
of ozone gain should be increased to 1 .2 to 2%. If 
the details of the stratospheric chemistry outlined 
in this report are correct, we predict substantial 
ozone gains from NOx emissions in the stratosphere. 
This prediction differs radically from those reported 
before 1977, all of which suggested that substantial 
ozone losses would be caused by the operation of 
supersonic aircraft, but ozone increases would occur 
when subsonic emissions were considered. Con- 
versely, water-vapor emission, which was once 
believed to have negligible effect on stratospheric 
ozone, may cause the destruction of a substantial 
amount of ozone with no NOx emission; perhaps an 
ozone loss of as much as 1% could result from a 
fleet of 250 SST’s, although the true value would 
probably be less, depending on the magnitude of the 
thermal feedback effect. A combination of NOx 
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H 2 O emission results in a small net ozone increase. 
Finally, and most important in the current context, 
the predicted latitudinal variation of the ozone 
changes using the Ames model has undergone marked 
modification. Determination of the meridional 
variation of ozone change is important because the 
amount of solar UBV that penetrates the ozone 
shield is dependent on the local optical depth of the 
ozone on the slant path to the Sun. We no longer 
predict the largest effect in the Northern Hemisphere. 
For aircraft operating at a nominal cruise altitude 
of 20 km, we find only a variation between equator 
and pole of ^^50% if water vapor is not effective in 
destroying ozone (fig. 29), even though most of the 
NOx is injected at latitudes between 30° and 60° N. 
The small interhemispherical asymmetries shown in 
figure 29 are due to differences in transport. If water 
vapor is effective, we find that ozone gains may be 
smaller in the Northern Hemisphere than in the 
Southern Hemisphere. However, if the operating 
altitude is lowered to 17.5 km, the meridional varia- 
tion is very small. Of course, for an SST fleet of 
about 250 aircraft, the ozone gain due to existing 
fleets of subsonic aircraft would far outweigh the 
supersonic aircraft effects (cf. fig. 24, 25). Inclusion 
of washout-rainout effects may also alter the results 
of scenarios (4) and (5). 

Despite the recent advances in stratospheric 
aeronomy discussed in this report, some significant 


uncertainties remain. Because of the predictions of 
HNO 3 /NO 2 ratios that are apparently too large 
(compared with measurements — see fig. 11 ), we may 
have overestimated the low-altitude ozone increase 
due to reaction ( 1 ); it is thus conceivable that air- 
craft operations in the stratosphere could still lead to 
small ozone losses. The importance of thermal 
feedback from the water-vapor injection into the 
ozone chemistry is very uncertain. This thermal 
coupling, which is due to infrared cooling by the 
water vapor deposited in the stratosphere, has the 
effect of slowing some of the reactions that destroy 
ozone. As we have emphasized several times in our 
discussion, large uncertainties in the transport rates 
of trace substances in the stratosphere cause con- 
siderable uncertainty in computed ozone change 
because of the cancellation effect of low-altitude 
ozone production and high-altitude ozone destruc- 
tion. Finally, the effects of ozone redistribution 
(increase below ^^25 km and decrease above) 
especially on stratospheric thermal structure, have 
never been assessed. All of these uncertainties apply 
not only to the globally averaged ozone change, 
but to meridional variations as well. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, California 94035, 

August 22, 1980 
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TABLE I 1990 WORLDWIDE AIRCRAFT NOx EMISSIONS, 
HIGH ESTIMATES, kg/yr" 


Latitude 




NOx eniissions, kg/yr 




6-8 km 

8-9 km 

9-lOkm 

10-11 km 

lM2km 

12-13 km 

13-14 km 

14-15 km 

N 60+ 

3.35E6 

3.03E6 

1.43E7 

131E7 

1.46E7 

1.31E6 

9.99E5 

2.03E5 

50-60 

2.15E7 

2.59E7 

9.44E7 

1.06E8 

9.09E7 

8.26E6 

3.72E6 

1 .06E6 

40-50 

7.60E7 

8.70E7 

1.79E8 

2.79E8 

1.62E8 

2.48E7 

4.59E6 

1.04E6 

3040 

7.74E7 

9.20E7 

1 .67E8 

3.09E8 

1 .72E8 

2.97E7 

3.1 1E6 

8.70E5 

20-30 

2.61E7 

2.83E7 

6.74E7 

1 .02E8 

6.92E7 

8.73E6 

1.55E6 

6.00E5 

10-20 

1.11E7 

1.18E7 

2.65E7 

4.28E7 

3.99E7 

3.67E6 

4.74E5 

7.70E4 

0-10 

4.80E6 

5.14E6 

1 .50E7 

1 .82E7 

1 .36E7 

1.26E6 

1 .73E5 

0 

10-0 

3.31E6 

3.77E6 

1 .22E7 

1 .38E7 

1 .09E7 

8.65E5 

1 .38E5 

0 

20-10 

2.74E6 

3.21E6 

U4E7 

1 .52E7 

1.15E7 

1.1 1E6 

3.15E5 

6.60E4 

30-20 

3.67E6 

4.01 E6 

9.47E6 

1.37E7 

8.66E6 

9.31E5 

5.10E4 

0 

40-30 

4.01 E6 

4.63E6 

6.62E6 

1.1 8E7 

6.14E6 

1.21E6 

8.64E4 

2.58E4 

5040 

2.36E5 

3.05E5 

3.19E5 

8.28E5 

4.46E5 

9.29E4 

1.5 El 

0 

60-50 

4.77E4 

3.79E4 

2.99E4 

2.52E4 

1 .04E4 

1.45E3 

0.97 

0 

S 60+ 

0 

0 

0 

0 

0 

0 

0 

0 

Total 

2.343E8 

2.69 1E8 

6.036E8 

9.255E8 

5.999E8 

8.197E7 

1.521E7 

3.94E6 


^Oliver et al. (1977) 
*2J43E8= 2.343X10® 
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TABLE 2.- COMPARISON OF THE AMES AND AEROSPACE TWO-DIMENSIONAL MODELS 


Characteristic 

Ames model 

Aerospace model 

Vertical dimensions 

0 to 60 km 

0 to 50 km 

Horizontal dimensions 

80° S to 80° N 

90° S to 90° N 

Vertical resolution 

2.5 km 

1 km from 0 to 35 km, 2.5 km from 
35 to 50 km 

Horizontal resolution 

5° 

10 ° 

Transport, advective 

Parameterized 

Parameterized 

Transport, eddy 

Parameterized 

Parameterized 

Transport, temporal changes 

Seasonal 

Monthly 

Upper boundary conditions 

Diffusive equilibrium 

Photochemical equilibrium 

(0,0CD),03,0H,H02,N,H) 
flux condition (other species) 

Lower boundary conditions 

Chemical equilibrium (0, OH, 
HO 2 , and other reactive 
species) 

Chemical equilibrium (0, OH, HO 2 , 
and other reactive species) 

Includes rainout-washout 

No 

Yes (for NOx, H 2 0, H 2 O 2 , HNO 3 ) 

Includes methane oxidation 

Yes 

Yes 

sequences 



Includes chlorine 

, Yes 

No 

Transports water vapor 

Yes 

Yes 

Reaction rates 

i See tables 3-7 

See tables 3-7 

Diurnal averaging 



Photo rates 

Yes 

Yes 

Rate coefficients 

Yes 

Yes 

Temperature feedback 

No 

No 

Time step 

1 day 

2 day 
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TABLE 3 .- REACTION-RATE CONSTANTS FOR THE 0 -N-H SYSTEM 


Reaction 

Rate constant^ 

Note 

Reaction 

Rate constant^ 

Note 

O-i-Oj+M-^Oj-tM 

1.1X10"®'^ 

b,c 

N + O3 ^ NO + O2 

2 . 0 X 10 “** 6“**“°/^ 

c 


1 . 07 X 10 '^'* 



5 . 0 X 10 “*^ 


O + O + M-^O^-tM 

3. OX 10 ’^^ ( 300 /T)^ 

b 

N + NO->0-tN2 

8.2X10“** 6“"**°''^ 


O + Os^Oj +02(‘SJ) 
0 + OH ^ H + O2 

1.9X10-‘‘ 

b 

N + O2 ^ NO + 0 

5.5X10"*^ e- 322 o/r 


4 . 2 X 10 ““ 


N+NO2 ^NjO-tO 

2.0X10“** e'®**®/^ 

d 

O3 + H OH + O2 

1.2X10“*° 


NO2 +03-^ NO3 + O2 

1.2X10“*® 


0 + HO2 ^ OH O2 

3 . 5 X 10 “** 


NO2 + 0 + M NO3 H- M 

1.0X10"®* 

b, e 

O3 + OH -> HO2 + O2 

1 . 5 X 10 “*^ 


NO + NO3 -> NO2 + NO2 

8 . 7 X 10 “*® 

b 

H + O2 + M -> HO2 + M 

2.1 X 10 “°^ 


NO2 +NO3 +M^N20s +M 


by e 

NO -t O3 ^ NO2 + O2 

2.1X10“*^ e-i 45 o/r 


N2O5 -t M ->■ NO2 + NO3 -t M 


b, e 

O + NO2 ^NO-t-02 

9.1X10“*^ 


NO + HO2 ^ NO2 + OH 

i 8 . 0 X 10 “*® 


NO + 0 + M NO2 + M 

1.6X10“°^ e=»^/^ 


NO2 + OH -t M ^ HNO3 + M 


f 

OH + OH->H 20 + 0 

1.0X10'** 


HNO3 -t 0 ^ OH -t NO3 

1 . 0 X 10 “*“* 

b,g 

OH-tHOj ->H2 0 + 02 

3 . 0 X 10 “** 


HNO3 +0H^H20 + N03 

8.0X10“*“* 


HO2 + HO2 H2 O2 + O2 

2 . 5 X 10 “*^ 



1 . 8 X 10 “®® e**®®/^ 


H2O2 +0H-^H20 + H02 

1.0X10“** e“’ 5 °/^ 


NO -t OH + M ^ HNO2 + M 

b,g 

H2O2 + 0 ^ 0 H + H02 

2.8X10“*^ 

c 

HNO2 + 0 ^ OH + NO2 

1 . 0 X 10 “*“* 

b,g 


j 2 . 75 X 10 “*^ 


HNO2 +0H->H2 0 + N02 

8.0X10“*“* 

1 

O3 + HO2 ^ OH -t O2 + O2 

: 7 . 3 X 10 “*"* e“*^'**^^ 





0H + 0H + M^H2 02 -tM 

1 . 26 X 10 “°^ 

1 . 3 X 10 “°^ e^oolT 

c 



! 


'^Rate constants are in molecule •cm* sec units; unless otherwise noted, the rate constants are taken from Hudson ( 1977 ). ^For 
detailed data reference see the reaction tabulation in Turco and Whitten ( 1975 ). ^The upper expression was used by Ames 
investigators, the lower one by Aerospace Corporation workers. ^This reaction was included in the Aerospace model only. 
^For the complete pressure-dependent rate expression, see Turco and Whitten ( 1975 ). •'The expression adopted for this rate 
constant is given in Hudson ( 1977 ). ^These reactions were included in the Ames model only. 



TABLE 4.- REACTION-RATE CONSTANTS FOR 
EXCITED OXYGEN SPECIES 


Reaction 

Rate constant^ 

Note 

0CD)+N2 ->0 + N2 

2.0X10-“ 6“’/^ 

b 

O0D) + O2^O + O2('sJ) 

2.9 X 10'“ 

b 

0(‘D) + H20^0H + 0H 

2.3X10"“’ 


0('*D) + N2 0-»-N0 + N0 

5.5X10"“ 


0(‘D) + N20-»-N2 +O 2 

5.5X10"“ 


0('D) + H2 ^OH + H 

9.9X10"“ 


0CD) + N2 +M->-N20 + M 

3.5X10"“ 



^Rate constants are in mo]ecule-cm*sec units. Unless otherwise 
noted, the rate constants are taken from Hudson (1977). 

^These reactions were included in the Ames model only. 



TABLE 5.- REACTION-RATE CONSTANTS FOR 
CARBON COMPOUNDS 


Reaction 

Rate constant^ 

Note 

CO + OH -> CO2 + H 

1.4X10"1^ 

b 

CH4 +0H^CH3 +H2O 

2.4X10-12 e-i’2 0 /r 



2.36X10-12 e- 1210 /r 


CH4 + 0 ^ CH3 + OH 

2.8X10-11 

c 


3.5X10-11 


CH4 + 0 (’ D) ^ CH3 + OH 

1.3X10-1° 


CH4 + 0 (*D)^CH 2 0 + H2 

1.4X10-11 

b 

CH3 +O2 +M->CH302 +M 


dy e 

CH3 +0->CH20 + H 

1.2X10-1° 

d, e 

CH 20 + 0 ->CH 0 + 0 H 

2.0X10-11 e-‘^®°/^ 


CH20 + 0H-^CH0 + H20 

3.0X10-11 e"25°/^ 



1.4X10-11 


CHO + O2 ^ CO + HO2 

6.0X 10-12 


CHO + 0 ^ CO + OH 

1.0X10-1° 

d, e 

CHO + 0 ^ CO2 + H 

7.3X10-11 

d,e 

CH3O2 +N0^CH3 0 + N02 

2.0X10-12 

bj 


1.5X10-12 


CH3O2 +CH3O2 ->2CH3 0 + 02 

2.0X10-1® 

d, e 

CH3O + O2 ->CH2 0 + H02 

1.6X 10-12 g- 3300 / 7 ’ 


CH3O2 +HO2 ^CH4 02 +O2 

2.5X10-12 

d 

CH4O2 + 0 H^CH 302 +H2O 

l.OX 10-11 g- 750 /T’ 

d, e 


^Rate constants are in molecule*cm*sec units; unless otherwise 
noted, the rate constants are taken from Hudson (1977). ^The 
upper expression was used by Ames investigators, the lower 
one by Aerospace Corporation workers. ^These reactions were 
included in the Aerospace model only. ^For detailed data refer- 
ences, see the reaction tabulation in Turco and Whitten (1975). 
^These reactions were used in the Ames model only.-^Based on a 
lower limit measurement of Cox (1977). 




TABLE 6.- PHOTODISSOCIATION RATES FOR 0-N-H-C 
CONSTITUENTS 


Photodissociation process^ 

Dissociation rate,^ 
sec”^ 

Oj + ht;— — ^ D) + 0 

6 . 6 X 10"’ 

^ . . 176<X<242nm „ . ^ 
O 2 + hp — ■ — ■■■»>0 + 0 

5.9X10"® 

O 3 + nm ^ 

2.1 X lO"* 

O 3 + h^^<310nin ^ 

3.9X10"® 


8.1 XIO"'^ 

^ , 450<X<750nm ^ ^ 

O 3 + ht^ — ■ ' »Oa + 0 

2.3X 10"^ 

NO + hi; ^ N + 0 

6 . 1 X 10 '® 

NO 2 + hi^ NO + 0 

6.4X 10'® 

NO 3 + hi/ “H. NO + 0 

2 . 0 X 10 "’ 

NO 2 + 0 

5.0X10"’ 

N20 + hi/^N 2 +0(‘D) 

9 . 8 X 10 "’ 

N 2 O 5 +hi/^N 02 +NO 2 +0 

3.1X10"® 

HNO 2 + hi; ^ OH + NO 

4.8X 10"® 

HNO 3 + hi; OH + NO 2 

9.0X 10'® 

H 20 + hi;^ 0 H + H 

4.4X 10'® 

H 2 O 2 + hi; ^ OH + OH 

7.1X10'® 

HO 2 +hi;-^0H + 0 

4.4 X 10"® 


6.4X10'®,(c) 

CH 20 + hz;^CH 0 + H 

6 . 6 X 10"® 

CH 2 0 + hi; -► CO + H 2 

8.9X10"® 

CH 4 O 2 +hi;->CH 30 + 0 H 

7.1X10"® 


^For detailed data references and discussion, see Turco (1975) 
and Turco and Whitten (1975). 

^Twenty -four hour average photodissociation rates at 120 km 
altitude are given. 

^This process was included in the Aerospace model only. 




TABLE 7.- PHOTOCHEMICAL REACTIONS AND RATE 
COEFFICIENTS FOR CHLORINE COMPOUNDS 


Reaction 

Rate constant^ 

Note 

Cl + O 3 ^ CIO + O 2 

2.7X10"" 


Cl + CH 4 ^ HCl + CH 3 

7.3X10"" 


Cl + H 2 ^ HCl + H 

3.5X10"" 


Cl + HO 2 ^ HCl + O 2 

3.0X10"" 


CIO + 0 ^ Cl + O 2 

7.7X10"" 


CIO + NO Cl + NO 2 

1 . 0 X 10 ”" 


CIO + NO 2 + M ^ CIONO 2 + M 


b 

CIONO 2 + 0 ^ CIO + NO 3 

3.0X10'" 


Cl + CIO 2 ^ CIO + CIO 

5.9X10"" 


HCl + OH^Cl + HzO 

3.0X10"" 


HCl + 0 -> Cl + OH 

1.1X10'" 


C10 + hi^->Cl + 0 

2 . 8 X 10 '^ 

c 

CIONO 2 +hi^^CI0 + N 02 

5.6X10''* 

c 

HCl + -> Cl + H 

1 . 6 X 10 '® 

c 

CF 2 CI 2 +hi^->Cl + Cl 

1.3X 10 '® 

c 

CFCI 3 +hi>^Cl + Cl + Cl 

7.1X10”® 

c 


^Rate constants are in molecule 'em* sec units; unless other- 
wise noted, the rate constants are taken from Hudson (1977). 

^The expression adopted for this rate constant is given in 
Hudson (1977) as computed by Ames workers. 

^Twenty -four hour average photodissociation rates at 
120 km altitude are given. 





40S 20 0 20 40 60 SON 

LATITUDE, deg 


Figure L - Nitrogen oxide emission rate in molecules of NO cm"^ sec'^ which represents the SST emissions model 
discussed in the text; the mismatches among the three curves occur because of the complex nature of the 
modeled air route structure and the increase in fli^t altitudes as fuel is consumed. 
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INTERPRETATION OF INITIAL 

DATA BY JOHNSTON et al. (1976) 



Figure 2. — Comparison of calculated and observed carbon-14 profiles at lat. 30'" N 9 months after an atmospheric 
nuclear test explosion; observational data are from Telegadas (1971). 
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AEROSPACE O SCHMELTEKOPF et al., (1977) 

ARC □ EHHALT (1978) 

A VEDDER etal., (1978) 



0° 40° 70° 



10-9 1Q-® 10-’ 10"® (0°) 


MIXING FRACTION 


Figure 3. — Comparison of calculated and measured profiles of nitrous oxide (N 2 O) at three latitudes in the Northern 

Hemisphere, (a) Summer, (b) Autumn. 
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Figure 4 . - Comparison of calculated and measured profiles of methane (CH4) at three latitudes in the Northern 

Hemisphere, (a) Summer, (b) Autumn. 
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io>^ ^ - — r z u \ I V 

80S 60 40 20 0 20 40 60 SON 


LATITUDE, deg 
(a) Summer. 

Figure 5.- Computed isopleths of ozone concentration in units of 10^^ cm"^ ; the solid lines represent the Ames 

model, the broken lines the Aerospace model. 
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AEROSPACE 

ARC 




10 ^^ 10^2 ( 0 °) 

CONCENTRATION, cm"^ 

Figure 6.- Comparison of calculated and measured profiles of ozone at three latitudes in the Northern Hemisphere. 
The squares represent the mean seasonal measured values obtained by Wilcox et al. (1977). The Goddard data 
(triangles) (Private communication from I. Eberstein, Goddard Space Flight Center [1977]), were obtained 
from satellite measurements of backscattered solar ultraviolet radiation (Heath et al. 1973). (a) Summer, 
(b) Autumn. 
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Figure 7.— Computed (using the new value for the rate coefficient of reaction (7)) and measured ozone concentra- 
tion profiles for autumn in the Northern Hemisphere. The squares represent the mean seasonal measured values 
obtained by Wilcox et al. (1977). The Goddard data (triangles)(Private communication from I. Eberstein, 
Goddard Space Flight Center [1977]) were obtained from satellite measurements of backscattered solar ultra- 
violet radiation (Heath et al., 1973). The solid lines correspond to the older value for the rate coefficient of 
reaction (9) and the broken lines to the value reported by Zahniser and Howard (1978). The Ames model (broken 
line) simulated 2.5 ppbv Clx; the Aerospace model simulated no Clx, (a) Equator, (b) Latitude 40° N. 
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LATITUDE, deg 


Figure 8. - Total ozone vertical column densities in the Northern Hemisphere. The solid lines represent computed 
values corresponding to the older value of the rate coefficient of reaction (9). The broken lines represent the 
computed values that correspond to the value reported by Zahniser and Howard (1978). The squares repre- 
sent the mean seasonal values obtained by Wilcox et al. (1977). (a) Summer, (b) Autumn. 
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NOy -SUMMER (NH) 



10-11 1 Q -10 10-9 10-8 ^ 0 ®) 

MIXING FRACTION 


NOy -AUTUMN (NH) 



Figure 9. — Computed profiles of odd-nitrogen, NOy, mixing ratios at three latitudes in the Northern Hemisphere. 
The data points are from a compilation by Ackerman (1976) of various measurements of NO, NO 2 , and HNO 3 
mostly at mid-latitudes; they are shown in figure 9 (a) for convenience, (a) Summer, (b) Autumn. 
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(a) Ames model. 

Figure 10.— Nitric acid column densities above 12.5 km as a function of latitude and season; balloon measurements 

were made in 1970. 
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HNO 3 /NO 2 AMBIENT CONCENTRATION RATIOS 
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Figure 11 . — Comparison of measured and computed values of the HNO 3 /NO 2 abundance ratio; the Harries (1978) 
data are from a re-analysis of the measurements first reported by Harries et al. (1976). 
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Figure 12. - Computed profiles of water-vapor volume mixing ratios at three latitudes in the Northern Hemisphere. 
The mean mid-latitude Northern Hemisphere vertical distribution of water vapor derived by Harries (1976) 
from available data are also shown; the dotted lines represent the rms variance among the individual measure- 
ments at each level. The measured data are plotted in figure 12 (a) for convenience, but they include measure- 
ments made in all seasons, (a) Summer, (b) Autumn. 
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Figure 13. — Comparison of calculated and measured vertical distributions of the hydroxyl radical OH; the measured 
values were obtained by balloon flights at Palestine Texas, lat. 32°(N.) (a) Summer (measured in July 1975; 
Anderson, 1976). (b) Autumn (measured in January 1976; Anderson, 1976). 
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Figure 14. — Comparison of calculated vertical distributions of the hydroperoxyl radical HO 2 . (a) Summer. 

(b) Autumn. 
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Figure 16. - Comparison of calculated and measured vertical distributions of NO^ (= NO + NO 2 ). The measured 
values were obtained over Aire sur I’Adour, France (lat. 43°35'N) in May 1974 (Ackerman et al., 1975). 

(a) Summer, (b) Autumn. 





Figure 17. - Odd-nitrogen increments (ANOy) that correspond to scenario (1). (a) Summer, (b) Autumn. 





Figure 18. - Odd-nitrogen increments (ANOy) that correspond to scenario (2). (a) Summer, (b) Autumn. 
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Figure 19.— Concluded. 





Figure 20. — Odd-nitrogen increments (ANOy) that correspond to scenario (4). (a) Summer, (b) Autumn. 
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Figure 22.— Absolute ozone changes (AO 3 ) that correspond to scenario (1) 
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Figure 23 - Absolute ozone changes (AO 3 ) that correspond to scenaiio (2), 







5] 


40 


AUTUMN 

ARC 

AEROSPACE 

INJECTION AT 20 km 

AND SUBSONIC AIRCRAFT 
5 years OF OPERATION 



10 


OL- 

-4 



(b) Au 


Figure 24.* 
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Figure 25 — Concluded. 
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Figure 27 — Absolute ozone changes (AO 3 ) that correspond to scenario (4), but without water-vapor injection 
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Figure 27.- Concluded. 
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Figure 28. 
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Figure 29.— Predicted ozone column change as a function of latitude for scenario (4), but without water-vapor 

injection. 
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Figure 31. 
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Predicted ozone column change as a function of latitutde for scenario (4) with the new value of the 

rate coefficient for reaction (2). 
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Figure 32. — Predicted ozone column changes as a function of latitude for scenario (5). 
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